Ballistic Spin Injection from Fe into ZnSe and GaAs 
with a (001), (111), and (110) orientation 
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We present first-principles calculations of ballistic spin injection in Fe/GaAs and Fe/ZnSe junc- 
tions with orientation (001), (111), and (110). We find that the symmetry mismatch of the Fe 
minority-spin states with the semiconductor conduction states can lead to extremely high spin po- 
larization of the current through the (001) interface for hot and thermal injection processes. Such a 
symmetry mismatch does not exist for the (111) and (110) interfaces, where smaller spin injection 
efflciencies are found. The presence of interface states is found to lower the current spin polarization, 
both with and without a Schottky barrier. Finally, a higher bias can also affect the spin injection 
efficiency. 

PACS numbers: 72.25.Hg, 72.25.Mk, 73. 23. Ad 



I. INTRODUCTION 

The idea of exploiting both the charge and the spin of 
an electron in semiconductor (SC) devices has lead to the 
growing field of spintronicsi Several spintronic devices 
are already proposedjSi^ but the full potential of spin- 
tronics has still to be discovered. At the moment the es- 
sential ingredients for spintronics, e.g., the injection and 
detection of a spin polarized current, the spin transport 
in SC etc., are only partially achieved in experiments. 
In this paper we want to investigate one of these main 
challenges: How to create a spin polarized current in a 
SC? There are several methods that obtain a high spin 
polarization, but most are useful for basic studies rather 
than applications: illumination of the SC with circularly 
polarized light, '^''^ injection of electrons from a ferromag- 
netic scanning tunneling microscope (STM) tip6'7 or from 
a paramagnetic semiconductor polarized in an external 
magnetic fieldAS Better applicable for devices is the in- 
jection from a diluted magnetic semiconductor^2*ii but 
these have the drawback of low Curie temperatures Tc up 
to now. More promising is the injection from a common 
ferromagnetic metal, like Fe, into a SC due to the high 
Tc- In this case it has been showniS. that in the diffusive 
limit the spin polarization of the injected current is van- 
ishingly smallf2iiii4 if the contact between the FM and 
the SC is ohmic. This 'fundamental obstacle' is traced 
back to the conductivity mismatch between both mate- 
rials. It can be overcome by either a ferromagnet with 
nearly 100% spin polarization at the Fermi energy, e.g., a 
half metal, or by inserting a spin polarized interface resis- 
tance at the FM/SC interfacefi^*i& e.g., an extrinsic tun- 
neling barrier or an intrinsic Schottky barrier. By means 
of a Schottky barrieriiLiSiiS or of an AI2O3 barrier— spin 
injection could be achieved in experiments. 

Here we investigate the spin injection from a ferro- 
magnet (FM) into a SC in the ballistic regime, i.e., the 
limit where inelastic and incoherent scattering events are 
negligible. In this case there is no scattering in the per- 
fect bulk regions and the whole resistance results only 
from the interface region so that there is no conductivity 



mismatch. Therefore other effects emerge that are not 
present in the diffusive limit: KirczenowSi has shown 
that some FM/SC interfaces can act as ideal spin fil- 
ter. This is the case, if for a direct-gap SC the Fermi 
surface of the FM projected onto the two-dimensional 
Brillouin zone (2DBZ) has a hole at the F point for one 
spin direction. But the material combinations Fe/ZnSe 
and Fe/GaAs investigated here, do not fulfill this condi- 
tion. It has been pointed out that the symmetry of the 
FM wave functions is very important for the spin injec- 
tion process—iSiS^ If there is a second FM lead to detect 
the spin polarization of the current, additional effects like 
the forming of quantum well states, Fabry-Perot like reso- 
nances and extremely high magnetoresistance ratios have 
been reported in a recent theoretical ab-initio study— 

There are two types of calculations for the ballistic 
spin injection already published: by means of analyt- 
ical models using plane wavesS^s^SiSi and by ab-initio 
methods— iSSiSiS^ While the model calculations obtain 
only a few percent of spin polarization, the inclusion of 
the full band structure of the FM and the interface region 
in ab-initio methods can result in spin polarizations up 
to 99%. The origin of this high polarization are the dif- 
ferent symmetries of the Fe d states at the Fermi energy 
for the majority and the minority spin. 

Here we extend our previous workSSiS^ and show more 
details of the spin injection process. In Sec. |n]the details 
of the calculations and the investigated Fe/SC systems 
are given. In Sec. IIIII the ballistic spin injection of hot 
electrons for the Fe/SC(001), (111), and (110) interfaces 
is investigated. It emerges that the (111) and (110) ori- 
entations do not show the large symmetry enforced spin 
polarization as the (001) orientation does. In Sec. llVl we 
report on injection of thermal electrons with and without 
a Schottky barrier, addressing also the effect of resonant 
interface states. Sec. gives an approximation to the 
effect of a higher bias voltage (out of the linear-response 
regime) . The paper closes with a summary of the results 
in Sec. EH 
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II. CALCULATION DETAILS 



The investigated heterostructures consist of a half- 
space of ideal bulk Fe and one of the ideal bulk SC (ei- 
ther ZnSe or GaAs). Their properties are described by 
the surface Green's function determined by the decima- 
tion technique^ Between the half-spaces there are sev- 
eral monolayers (MLs) of the interface region, where the 
atomic potentials are allowed to deviate from their bulk 
values. This interface region consists of four MLs of Fe 
and two MLs of SC for the (001) orientation and appro- 
priate numbers of MLs for the (111) and (110) orienta- 
tions. The corresponding potentials are obtained self- 
consistently from a Fe/SC/Fe junction geometry by the 
screened Korringa-Kohn-Rostoker (KKR) Green func- 
tion methodi^fliSi It is assumed that the SC lattice is 
matched to the bulk Fe, resulting in an fee SC lattice con- 
stant double the lattice constant of bcc Fe dsc = Sdpc = 
5.742A. The experimental lattice constant of Fe is used. 
This means that the SC is stretched compared to the bulk 
lattice constant by 1.3% and 1.6% for ZnSe and GaAs, 
respectively, resulting in a slightly smaller energy gap. In 
order to describe the SC accurately, two empty spheres 
per unit cell are introduced to account for the open space 
in the geometry of the zinc-blende structure. A perfect 
two-dimensional translation symmetry in the whole sys- 
tem is assumed and so the in-plane component k|| of the 
k vector is conserved while crossing the interface. The z 
axis is always assumed to be the growth direction, i.e., 
standing perpendicular on the interface. Also the con- 
ductance G is calculated by a two-dimensional Fourier 
transformation in dependence of k|| . 

The calculations of the ground state properties are per- 
formed within the density-functional theory in the local 
density approximation (LDA) for the exchange and cor- 
relation terms. It is well known that it underestimates 
the energy gap in SC by around one half. Otherwise it 
gives very accurate FM and SC bands and is well suited 
for the considered problem. The potentials are described 
within the atomic-sphere approximation, and the wave 
functions are expanded in angular momentum up to a 
cutoff of £max = 2 for self-consistency and imax = 3 for 
the conductance calculations. Spin-orbit coupling and 
spin- flip scattering are not included in the calculations. 
This could have an effect in case of injection into the SC 
valence band, where the spin-orbit interaction is stronger, 
but not in the case of the valence band which has mainly 
s character and negligible spin-orbit coupling. 

The conductance G is calculated by the Landauer 
formula. G is expressed as a sum of the transmis- 
sion probabilities through the scattering region (here the 
Fe/SC interface) over all available conducting channels. 
These channels are enumerated by the spin orientation 
a —1 or I, the bands v,!/' at the Fermi energy of the 
left and right lead, respectively, and the k|| vector. The 
spin-dependent conductance, normalized to the area of 



the two-dimensional unit cell, reads 

h A2DBZ , J 

2DBZ 

with A2DBZ the area of the two-dimensional Brillouin 
zone. Since the ky vector and the spin a of the electron 
are conserved in our calculation, they are the same for 
the incoming and the outgoing states. Here the Landauer 
formula is evaluated by a Green's function formalism in- 
troduced by Baranger and Stonci^ Details are presented 
elsewherei^ The spin polarization P of the current is de- 
fined by the spin dependent conductances in both spin 
bands 

G^_G^ 

GT + ' 

where G^ and G^ is the conductance of the majority and 
the minority electrons, respectively. 

III. HOT ELECTRON INJECTION PROCESS 

To investigate the effects of the symmetry of the Fe 
bands, we calculate first the hot injection of electrons, 
i.e., of Fe electrons well above the Fermi level Ep, which 
falls in the gap of the SC. For most of the calculations we 
restrict the discussion to the F (k|| = 0) point for sim- 
plicity, i.e., we consider only electrons with perpendicular 
incidence on the interface. This is motivated by the fact 
that in most spin injection experiments the Fermi energy 
in the SC is only some tens of meV above the conduction 
band minimum resulting in a very small Fermi sphere 
around the F point. 

A. The Fe/SC (001) orientation 

First we investigate the spin injection in Fe/SC junc- 
tions grown in the (001) orientation. The (001) direction 
has the highest symmetry and is therefore the best can- 
didate for a symmetry enforced high spin polarization of 
the injected current;^ i.e., almost total reflection of the 
minority spin electrons due to the symmetry mismatch 
between the SC and minority spin Fe bands. 

For the investigation of the bands that are available 
in Fe and in the SC at the F point, it is important to 
note that the in-plane common two-dimensional unit cell 
is double in area than the one of bulk bcc Fe (containing 
two inequivalent Fe atoms), and the common 2DBZ has 
half the area of the one of bcc Fe. Thus additional bands 
are backfolded into the 2DBZ and the number of bands 
at F increases. A description of the backfolding is given 
in the Appendix. 

The Fe band structure at the F point along fc^ (F-H, 
i.e. the A direction), including the backfolded bands, is 
shown in Fig.nfor both spins. This direction has in bulk 
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TABLE I: Symmetry properties of energy bands in Fe [001] 
{Cav symmetry group) and their local orbital analysis. The 
backfolded bands are indicated by (backf.). The last column 
shows the coupling to the SC Ai conduction band {C2v sym- 
metry group) based on symmetry compatibility relations and 
local orbital form (d orbitals are localized and expected to 
couple poorly; orbitals are extended into the SC and ex- 
pected to couple well). 



Rep. in Fc Orbital analysis 



Coupling to SC Ai 



Ai s;pz;d^2 good 

A2' dj:y poor 

A2 d^2_y2 none (orthogonal) 

A5 Px;Py;d^z;dy^ none (orthogonal) 

Di (backf.) s;px\dyz\d^2_y'2 — d^2 poor 

D2 (backf.) d^2_y2 +d^2 poor 

D3 (backf.) py+p^; d^y + d^^ poor 




FIG. 2: Band structure at the F point of the 2DBZ for ZnSe 
(001) (left panel) and GaAs (001) (right panel). The numbers 
denote the symmetry of the bands along the A directionj2& d 
is the SC lattice constant. 




FIG. 1: Band structure of Fe(OOl) at the F point of the 2DBZ 
(black lines represent the bulk F-H bands and gray lines the 
backfolded N-P-N bands). The left panel shows the majority 
and the right one the minority spin band structure. The Fermi 
energy is shown by a horizontal line. The numbers denote the 
symmetry of the bands along the A direction (F-H) and the 
D direction (N-P-N) 



Fe a symmetry, just as the Fe (001) sm'face. In the 
2DBZ of the SC no bands are backfolded to the F point. 
The ZnSe and GaAs band structmes along F-X (the A 
direction) are shown in Fig. |21 In the SC this direction 
has a C2V symmetry (a subgroup of C4„), as the SC (001) 
surface. Therefore the Fe/SC(001) interface has also the 
reduced C2V symmetry. As a result, some bands d that 
are mutually orthogonal in bulk Fe are able to couple to 
each other at the interface. 

Now it is important to examine which bands of Fe are 
allowed to couple to the Ai conduction band of the SC 
which transforms fully symmetrically under the rotations 
of C2v', this can be found by the compatibility relations 
between the groups €4^ and C2V Note that the sym- 
metry notation of Fe and the SC refers to the different 
groups, C4V for Fe and C2V for the SC. At Ep and above 
the available Fe bands along the A direction (black lines 
in Fig.^l and their symmetry and coupling properties are 



given in Tabled Practically, the only ones that are able 
to couple well are the Ai bands containing the extended 
s and Pz orbitals, as well as the (1^2 orbitals pointing into 
the SC, while the A2' bands couple only poorly because 
they consist of more localized d^y orbitals pointing in- 
planc. At Ep, and up to 1.2 eV above, the former bands 
exist only for majority spin, allowing transmission, while 
for minority spin the latter bands co-exist with the A2 
and A5 bands. 

In Fig. 13 the calculated conductance for the hot 
injection process for Fe/ZnSe(001) and in Fig. ^ for 
Fe/GaAs(001) are shown. In (a) the Zn and Ga termi- 
nated and in (c) the Se and As terminated interfaces are 
shown. The intermixed interfaces (b) will be discussed 
later. As it was already discussed above, the coupling of 
the minority A2' band is much weaker than the coupling 
of the majority Ai band to the Ai conduction band of 
the SC. Thus the A2' states are nearly totally reflected at 
the interface. This results for all three terminations in a 
very high symmetry enforced spin polarization. Around 
1.3 eV above the Fermi energy (at 11.3 eV) also a Ai 
band is available in the minority channel. At this energy 
the conductance in the minority channel rises drastically 
showing that the high spin polarization at lower energies 
originates precisely from the absence of the Ai Fe band 
in the minority channel. Also the flat backfolded bands 
of Di and D2 symmetries give a small contribution to 
the conductance in the small energy window of around 
0.3eV width, centered around 11.4 eV. At these energies 
in the ZnSe half-space there is only one conduction band 
available at the F point, limiting the maximum conduc- 
tance to 1 e^/h. In the case of the GaAs half-space there 
are three conduction bands and so the conductance could 
rise up to a maximal value of 3 /h. 

Finally the conductance can be seen to decrease for 
small conduction band energies and even vanishes at the 
conduction band minimum of the SC. This can be ex- 
plained qualitatively by the vanishing group velocity, so 
that no current can be transported away from the inter- 
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FIG. 3: Hot injection of electrons in Fe/ZnSe(001) with a Zn 
terminated (a) , an intermixed Zn terminated (b) and a Se ter- 
minated (c) interface. The black line shows the majority and 
the gray line the minority spin. The conductance is evaluated 
at the r point for simplicity. 
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FIG. 4: Caption as in Fig.Elbut for Fe/GaAs(001) with a Ga 
terminated (a), an intermixed Ga terminated (b) and an As 
terminated interface. Black lines refer to majority and gray 
lines to minority spin directions. 
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FIG. 5: Band structure of Fe(lll) at the F point of the 2DBZ 
(corresponding to the bulk band along the F-P-H symmetry 
line). The left panel shows the majority and the right panel 
the minority spin. The Fermi energy is shown by a horizontal 
line. The numbers denote the symmetry of the corresponding 
A- (between F and P) and F- (between P and H) symmetry 
directioni— 



metry with one atom per unit cell. But in the SC half- 
space this direction has only a three-fold rotational axis. 
In each ML there is either a cation or an anion (or one 
of the two types of empty spheres). Also there are a two 
kinds of possible geometric terminations: one where the 
terminated SC atom is singly and one where it is triply 
coordinated to the Fe atoms. 

In Fig.|Slthe band structure of Fe(lll) at the F point 
in the [111] direction is shown. In the SC(lll) 2DBZ no 
bands are backfolded to the F point and the only available 
band corresponds to the bulk band along the F-L high 
symmetry line. 

The band structure plot for Fe(lll) shows important 
differences compared to the (001) orientation: At the 
Fermi energy in both spin channels there are bands with 
the same Ai symmetry and there cannot be any symme- 
try enforced spin polarization. All eventually obtained 
differences in the conductances for both spin directions 
are due to the different coupling of the Fe bands to the 
SC. This coupling is believed to be sensitive to the exact 
interface properties like the termination, lattice relax- 
ations etc. The calculated results for the hot injection 
process in Fe/ZnSe(lll) with a Zn terminated interface 
are shown in Fig.El Also for this orientation the conduc- 
tance is shown at the F point as before. 

In both the singly and the triply-coordinated Zn ter- 
minations the obtained spin polarization is much smaller 
than in the (001) interface, demonstrating the lack of a 
symmetry enforced spin filtering. Above 12.0eV there is 
no ZnSe band exactly at the F point in the [111] direc- 
tion, and so no propagating states are available. 



B. The Fe/SC(lll) interface 

In this section the spin injection for the Fe/SC junc- 
tions grown in the (111) orientation is investigated. In 
the Fe half-space the [111] direction has a hexagonal sym- 



C. The Fe/SC(110) interface 

This section deals with Fe/SC(110) oriented junction. 
This is the one with the lowest symmetry considered in 
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Energy [eV] 

FIG. 6: Hot injection of electrons in Fe/ZnSe(lll) for a Zn 
terminated interface, where the Zn atoms are either triply 
coordinated (upper plot) or singly coordinated (lower plot) 
with respect to the interface Fe atoms. The black lines refer 
to the majority and the gray lines the minority spin direction. 
The conductance is evaluated at the F point for simplicity. 
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FIG. 8: Hot injection of electrons in Fe/ZnSe(110). The black 
hne show the majority and the gray line the minority spin 
conductance. Since all four different atoms (Zn, Se, and two 
vacancies) are located in each ML, there is only one possible 
termination. The conductance is evaluated at the F point for 
simplicity. 



the SC(llO) half-space there is also a backfolded band at 
the r point (created by applying bi in Fig. ll8|l . Because 
of this the maximum conductance per spin direction is 
2eV/i. 

The variety of backfolded bands and the reduced sym- 
metry of the [110] direction do not allow for a simple 
symmetry-based discussion. Basically, for the (110) ori- 
ented heterojunctions there is no symmetry enforced spin 
polarization. This is also shown in Fig.|Hl where the con- 
ductances at the r point for the hot injection process are 
presented. 

From the above we conclude that the best candidate 
for a symmetry-enforced high current spin polarization 
is the (001) interface. Thus, in what follows, we restrict 
our study to this case. 



FIG. 7: (color online) Band structure of Fe(llO) at the F 
point of the 2DBZ. The different colors denote the different 
backfolded bands corresponding to the following bulk bands: 
black: F-N, red: N-H, green: bi applied, blue: bi and b2 
applied (see the Appendix). The left plot shows the majority 
and the right plot the minority spin states. The Fermi en- 
ergy is shown by an horizontal line. The numbers give the 
symmetry of the corresponding bulk band;^ if possible. 



this work. From this point of view a high symmetry- 
enforced spin polarization is most improbable. Also the 
surface unit cell in the SC half-space contains all four 
kinds of atoms (cation, anion and both types of vacan- 
cies) and is larger than the (001) and the (111) unit cells. 
By matching the two dimensional translational symme- 
try, also the Fe surface unit cell has to contain four Fe 
atoms. This results in a very small 2DBZ for this ori- 
entation leading even in the SC half-space to backfolded 
bands at the F point. 

The band structure of Fe (001) at the F point with 
the different backfolded bands are shown in Fig. [7| In 



D. Hot injection process ky resolved 

In the hot injection process the restriction to the F 
point was motivated by the fact that in most applications 
the Fermi level is positioned only slightly in the conduc- 
tion band resulting in a very small Fermi sphere around 
the F point. More correctly the conductance should be 
integrated over the whole 2DBZ, because also states away 
from the F point are populated for higher injection ener- 
gies. Therefore in Fig. I^lthe k|| resolved conductances are 
shown for the majority and minority electrons for the Zn- 
terminated Fe/ZnSe (001) interface. The energy of the 
injected electrons is 500 meV above the conduction band 
minimum. In this case the polarization at the F point is 
99.7%. If the conductance is integrated over the 2DBZ 
the polarization is reduced to 39.5%. The reason is a 
higher conductance in the minority band for states away 
from the F point. The symmetry arguments discussed 
above are valid only at this high symmetry point. For 
states k|| ^ also other Fe minority states are allowed 
to couple to the SC Ai conduction band and so the con- 
ductance rises. We conclude that the high symmetry- 
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FIG. 9: Conductance in the hot injection process in Fe/ZnSe 
(001) (Zn termination), resolved in ky. The upper plot shows 
the conductance of the majority and the lower one for the 
minority spins. Only one tenth of the 2DBZ around the F 
point is shown. The Fermi energy lies 500 meV above the 
conduction band minimum in the semiconductor half-space. 
The reciprocal vectors are in units of 2n/da,c- 



enforced spin polarization can be realized only for ener- 
gies up to a few tens of meV above the conduction band 
edge. 

In the conductance plot (Fig. O for the minority elec- 
trons the reduced C2v symmetry of the Fe/SC(001) in- 
terface can be directly seen. The majority conductance 
behaves like free electrons traveling across a potential 
step and is practically unaffected by the reduced symme- 
try. The circularly symmetric form of 6*^(1411) and the 
twofold-symmetric form of G^(k||) can be understood in 
terms of the form of the Bloch functions for small devi- 
ations from k|| — in the same manner as described in 
Ref.Hfor the Fe/SC/Fe (001) junctions. 

E. The intermixed interface 

By means of an ab-initio study^^ it has been found that 
the interface of Fe/GaAs(001) with a Ga terminated in- 
terface is energetically more stable, if Fe atoms diffuses 
into the vacancy sites of the first Ga ML. In case of an 
As terminated interface the abrupt interface is found to 
be more stable than the intermixed one. The former case 



is accompanied by rather large lattice distortions which 
are not taken into account in our calculations. Although 
we do not know of analogous studies for Fe/ZnSc(001) 
interfaces we have also performed the intermixed inter- 
face in case of a Zn termination. The results are shown 
in Fig. inib) and 2{b) for the hot injection process with 
an intermixed Zn and Ga interface, respectively. The 
still high spin injection polarization can be explained by 
the fact that the intermixed Fe atoms do not change or 
reduce the C2V symmetry of the interface and so the sym- 
metry arguments are still valid. This is a different effect 
than interface roughness that reduces the spin polariza- 
tion drasticallyi^^i^ Solely the coupling of the Fe states 
to the SC are changed. In both cases the coupling in the 
minority band is enhanced in comparison to the atomi- 
cally abrupt interfaces shown in Fig.l^^a) and^a). 

IV. THERMAL INJECTION PROCESS 

In this section the injection of thermal electrons is in- 
vestigated. This is achieved by lowering the potentials 
in the SC half-space by a rigid shift, so that the Fermi 
energy falls slightly above the conduction band minimum 
(here around 10 meV). All other potentials are kept in 
their ground state position, in particular also the poten- 
tials of the two SC monolayers in the interface region. 
This shifting simulates the effect of n-doping or an ap- 
plied gate voltage in a field-effect-transistor device. The 
electrons are injected at the Ep of the Fe half-space di- 
rectly into the conduction band. Since the Fermi energy 
is only some tens of meV above the conduction band 
minimum and since effective masses of ZnSe and GaAs 
are small, the resulting Fermi wave vector in the SC half- 
space is very small (around one hundredth of the distance 
to the boundary of the 2DBZ). So we present here the 
conductance only for the F point. The evaluated po- 
larization changes only slightly, if the conductances are 
integrated over the whole 2DBZ. 



A. Thermal injection without Schottky barrier 

The results for Fe/ZnSe and Fe/GaAs(001) are dis- 
cussed in Ref. and show a very high spin polariza- 
tion of more than 97% for ZnSe and practically 100% for 
GaAs (001). Similar results are found for tunneling in 
the Fe/SC/Fe (001) junctionsi^^ Here we add the ther- 
mal injection with an intermixed Zn- and Ga-terminated 
interface (Fig. I10|l . In the case of an intermixed Zn in- 
terface nearly the same conductances are obtained as for 
the abrupt one (Fig. 3 in Ref. '23). But a drastic in- 
crease of the minority conductance can be seen for the 
intermixed Ga interface, resulting in a much smaller spin 
polarization. The reason for this increase is an interface 
resonance localized at the intermixed Fe atoms. This 
can be seen in the density-of-states (DOS) for this het- 
erojunction shown in Fig. II II for the intermixed Zn and 
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FIG. 10; Thermal injection of electrons in Fe/ZnSe (upper 
plot) and Fe/GaAs (001) (lower plot) with an intermixed Zn 
or Ga interface, respectively. The black line shows the major- 
ity and the gray line the minority spin. The conductance is 
evaluated at the F point. Note that the minority conductance 
is enlarged by a factor of 10 in case of the Zn termination. 
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FIG. 11: Density-of-states (DOS) for the thermal injection 
process in Fe/ZnSe (001) (upper plot) and Fe/GaAs (001) 
(lower plot) with an intermixed Zn or Ga interface. Arrows 
indicate the resonance at the intermixed Fe atom. 



Ga interface. The arrows indicate the position of the 
Fe resonance. In the case of the intermixed Zn interface 
this resonance lies slightly above the Ep and give rise 
to the small maximum in the minority conductance at 
Ep - Ec = 15 meV seen in Fig. UHl In the DOS for 
the intermixed Ga interface this resonance lies nearly on 
the Fermi energy and is much narrower, so that it has a 
strong influence on the minority conductance. Since the 
energy position of this resonance is believed to depend 
strongly on the interface properties, e.g., lattice relax- 
ations, further investigations are needed for the exact 
energy position. If the injection process is evaluated at 
only 0.3eV higher energies, the influence of this resonance 
is negligible and the normal high spin polarization is re- 



stored. However, for a higher bias one must integrate 
over an energy interval (see below) and thus include res- 
onant interface states that lie close to, but not at, Ep. 



B. Thermal injection through a Schottl4y barrier 



In 



theoretical 
as. 16 



studies for strongly diffusive 
transporWiiS it is shown that the negligible spin 
polarization due to the conductivity mismatch between 
the ferromagnetic metal and the semiconductor— can be 
overcome by a spin dependent tunneling barrier at the 
interface. Such a barrier could be the Scottky barrier 
create d at the interface. Although the approaches in 
Refs. Il2ll5ll()l assume diffusive transport, tunneling 
barriers can be described in the ballistic regime. 

In this section the thermal injection process is extended 
to the case where a Schottky barrier at the interface is 
included. The potentials for this junctions are the same 
as in the thermal injection process, but with an inserted 
tunneling barrier at the interface. The barrier is modeled 
by a rigid shift of each SC atom potential, constant within 
each monolayer. It starts at the third SC layer from 
the interface, and increases linearly in magnitude with 
distance from the interface. In this way the SC bands 
are gradually shifted downward. Finally, at the end of 
the barrier Ep lies 10 meV above the conduction band 
edge (and the shift is not changed from then on), while 
at the interface Ep lies in the middle of the gap. The 
barrier thickness is then varied between and 140 MLs. 

As shown in Ref. J2i a tunneling barrier at the Fe/SC 
(001) interface gives also a high spin polarization. Only 
in case of a resonant interface state near Ep, especially 
in Fe/ZnSe(001) with a Zn terminated interface, the po- 
larization is reduced. Here we will discuss the effect in 
more detail. In Fig. E| the conductance and the polar- 
ization in dependence of the barrier thickness is shown 
for Fe/ZnSe and Fe/GaAs(001), respectively. Due to the 
small Fermi wave vector, the conductance is again ana- 
lyzed only at the F point. The polarizations agree with 
the ones correctly integrated over the 2DBZ within 5%. 

Except for the Se terminated Fe/ZnSe(001) interface, 
all terminations show the influence of a resonant inter- 
face state in the minority band, resulting in an higher 
conductance for the minority spins with thicker barriers. 
This leads to a reduction of the spin polarization. For 
Fe/GaAs(001) this effect is also visible in the conduc- 
tance but is negligible for the spin polarization due to 
the much larger difference of the conductance for both 
spins than in Fe/ZnSe. To discuss this effect in more 
detail, in Fig. ^] the DOS at the F point is shown for 
Fe/ZnSe(001). The DOS for Fe/GaAs(001) is qualita- 
tively comparable for this purpose. In the DOS an inter- 
face state is visible near the Fermi energy in the minority 
spin channel. This state has Ai symmetry and lies in 
the energy region of the Ai hybridization gap of the Fe 
minority band structure. It can penetrate well a Schot- 
tky barrier of moderate thickness since the evanescent 



8 



Fe/ZnSe(001) 



Fe/GaAs(001) 

^ ■ ' ' • 

Ga term. 




50 100 

Barrier thickness [ML] 



50 100 150 

Barrier thickness [ML] 



FIG. 12: Influence of a smooth Schottky barrier on the 
spin dependent conductance (lower plots) and the spin po- 
larization (upper plots) for Fe/ZnSe(001) (left panel) and 
Fe/GaAs(001) (right panel). Circles show the conductance 
for the Zn or Ga and triangles for the Se or As termination. 
The conductance is evaluated at the F point. Arrows indi- 
cate the influence of the resonant interface state causing a 
non-exponential decay of the conductance for moderate bar- 
rier thickness. In Fe/ZnSe(001) the lines for the majority 
conductances lie on top of each other. 
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FIG. 14: Influence of the pinning position of the Fermi en- 
ergy relative to the middle of the gap on the conductance in 
Fe/ZnSe(001) junctions with a Zn terminated interface. The 
Schottky barrier is 80 ML thick. The black line shows the 
conductance for the majority and gray line of the minority 
spin. Negative energy means a pinning position of the Fermi 
energy near the valence band and positive ones near the con- 
duction band. The conductance is evaluated at the V point. 
The minority-conductance peak at —0.2 eV reflects the in- 
fluence of the interface state. 
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FIG. 13: DOS at the F point for Fe/ZnSe(001) with a Zn 
(upper plot) and a Se terminated interface (lower plot) . Gray 
lines show the DOS of the Fe interface ML and black lines the 
DOS of the Zn or Se interface ML. The vertical line indicate 
the position of the Fermi energy. All potentials are in the 
ground state position. Arrows indicate the peak in the DOS 
due to the resonant interface state. 



wave with the longest decay length is of the same Ai 
symmetryi^Si In the Fe half-space this interface state is 
normally also evanescent since there are no propagating 
states it can couple to. But due to the reduced symmetry 
of the interface the Ai interface state couples weakly to 
the Fe A2' band and becomes resonant. In the case of 
the Se terminated Fe/ZnSe(001) interface the interface 



state lies further away from Ep and practically does not 
contribute to the conductance. 

The potentials used for the DOS in Fig. E| arc the 
ground state potentials, i.e., Ep is in the middle of the 
gap in the SC half-space and no Schottky barrier is in- 
serted. If now a Schottky barrier is introduced into the 
junction, the system and the resonant interface state are 
slightly distorted, because the Ai interface state inter- 
acts with the Ai conduction band of the semiconductor. 
This results in a slight downshift in energy of the resonant 
interface state away from Ep for smaller Schottky barri- 
ers. The reduction of the spin polarization with thicker 
Schottky barriers can thus be explained by a shift of the 
resonant interface state towards Ep. 

This effect can also be seen by changing the pinning 
position of the Fermi energy relative in the SC gap at 
the interface (the potentials in the interface region are 
kept fixed). In Fig. ^] the conductance for different 
pinning positions in the gap is shown for Zn-terminated 
Fe/ZnSe(001) with a 80 ML thick Schottky barrier. The 
minority-conductance peak at —0.2 eV reflects the influ- 
ence of the interface state. 

Furthermore one can estimate the interface resistance 
under the assumption of diffusive transport in the bulk, 
using Schep's formula.--Mi Although this has been de- 
rived for metallic multilayers, we have applied it to this 
case. Our results are summarized in Table Ull More de- 
tails can be found in Ref. 
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P (Zn) 


P (Se) 


P (Ga) 


P (As) 


Ri [Qm^] 


(8 ML) 


1 • 10"^" 


8 • 10"" 


2 • 10~" 


2 • 10"^° 


P 


(8 ML) 


96% 


96% 


99.8% 


99.8% 


Ri [Qm^] 


(80 ML) 


5- IQ-^ 


7- 10"^ 


2 • 10"^ 


4 • 10"^ 


P 


(80 ML) 


77% 


97% 


99.2% 


98.0% 



TABLE IL Interface resistance Ri and spin polarization P 
for two Schottky barrier thicknesses and all four terminations 
in Fe/SC(001). Results are obtained by integration over the 
whole 2DBZ. A thickness of 80 ML corresponds to 115A. 



V. HIGHER BIAS VOLTAGE 

Up to now it was always assumed that the applied bias 
voltage is so small that only states at Ep carry current. 
But experimentally especially for thicker barriers, a bias 
voltage in the order of IV or less is used.^^ This relatively 
high bias in experiments with an optical detection of the 
spin polarization is needed to gain a reasonable signal-to- 
noise ratio of the emitted light. It is believed that the bias 
can be reduced in experiments with thinner tunneling 
barriers or with an electrical detection of the polarization, 
since then a much smaller current is needed. If a non-zero 
bias is applied, one must integrate the transmission over 
the appropriate energy range, and the Landauer formula 
reads^ 

-eU 


TT^E / d'hT^yih^Ep + E) (2) 

n A2DBZ , J 

"■'''■'^ 2DBZ 

with T^^, (k|| , Ep + E) being the transmission probability 
at the energy Ep + E. This equation gives the standard 
Landauer formula in the case of small bias voltages U 
in a first approximation, under the assumption that the 
electronic structure is unaffected by the applied voltage. 

Here we show that also for a non-zero bias and a tun- 
neling barrier of moderate thickness the T point is of main 
importance and a high spin polarization is obtained. It 
has been shown'^^ that in tunneling the smallest damp- 
ing factor K is located at the T point in the energy gap 
for ZnSe and GaAs (and for other direct-gap semicon- 
ductors), so that this point plays the major role for thick 
tunneling barriers. Also at this point the symmetry en- 
forced spin polarization is the highest, as shown in the 
k|| -resolved hot injection process. In the calculations the 
Fermi energy in the bulk SC is assumed to be 10 meV 
above the conduction band minimum. For the Zn termi- 
nation an abrupt interface is taken for simplicity. Due 
to the high bias voltage a k|| integration over the whole 
2DBZ is performed, because also states away from the F 
point in 2DBZ are populated and carry the current. 

The position of the energy bands and of the Fermi en- 
ergy in the Fe and SC half-spaces are sketched in Fig. 1151 




FIG. 15: Sketch of the position of the energy bands and the 
energies in the injection process with an applied bias voltage 
and a tunneling barrier. The bias voltage is IV. 



Barrier thickness (ML) 


P (Zn) 


P (Se) 


8 


63.6% 


61.4% 


80 


43.8% 


87.5% 



TABLE III: Current spin polarization in Fe/ZnSe(001) junc- 
tions with an applied bias voltage of IV. The results are ob- 
tained by an integration over the whole 2DBZ. In the Bulk 
ZnSe(OOl) the Fermi energy lies lOmeV in the conduction 
band. 



Already from this simple sketch it is qualitatively clear 
that the conductance at the upper energy range (at the 
Fermi energy of the Fe half-space) is the most important 
one, because there the effective tunneling barrier is the 
smallest. In Table Hill the spin polarizations are listed for 
a Fe/ZnSe(001) junction with an applied voltage of IV 
and a tunneling barrier of 8 ML and 80 ML thickness. 
In case of a Zn termination it is important that slightly 
below the Fermi energy of the Fe half-space the resonant 
interface state contributes largely to the minority conduc- 
tance. This gives a relatively low and a more or less con- 
stant spin polarization with thicker tunneling barriers. In 
the Se terminated interface the resonant interface state 
lies lower in energy, where the effective tunneling barrier 
is thicker than for the Zn terminated interface. Thus the 
resonant interface state has a much smaller influence on 
the minority conductance. With larger barrier thickness 
the high spin polarization is more or less restored because 
only states at F, having the smallest decay parameter are 
important and the symmetry arguments apply again. 



VI. SUMMARY AND CONCLUSIONS 

We have reported first-principles calculations on spin 
injection from Fe into GaAs and ZnSe through the (001), 
(111), and (110) interfaces. The electronic transport has 
been assumed to be in the ballistic regime, and the in- 
terfaces have been assumed to have two-dimensional pe- 
riodicity, so that k|| is conserved during scattering at the 
interface. Then, in a Landauer-Biittiker approach, the 
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conductance is determined by the transmission probabil- 
ity summed for all k|| . Under these assumptions, we have 
considered hot and thermal injection — the latter also in 
the presence of a Schottky barrier — and approached 
the higher bias regime. We have reached the following 
conclusions. 

(i) The spin polarization of the current is highest when 
most incoming Fe Bloch states of one spin direction are 
totally reflected due to symmetry mismatch with the SC 
conduction band states — we call this effect "symmetry- 
enforced spin polarization" . This is the case for the (001) 
interface, where minority-spin electrons for k|| = prac- 
tically cannot be transmitted. The lower symmetry of the 
(111) and (110) interface does not lead to such a selec- 
tion rule. Even for the (001) interface, in the case of hot 
injection, above an energy threshold where the Fe minor- 
ity Ai band starts, the minority-spin current increases 
rapidly even at k|| = 0, being then symmetry-allowed. 

(ii) For the symmetry-enforced spin polarization to be 
realised, the injection must take place close to the con- 
duction band edge, so that only states near k|| = are 
relevant; these have a well-defined and suitable symme- 
try properties. Also, the interface should be as ordered 
as possible, otherwise k|| is not conserved and incoming 
minority-spin Fe states from all k|| can scatter into the 
SC conduction band. Then the spin polarization will de- 
crease. 

(iii) The case of thermal injection (i.e. exactly at Ep), 
has been considered for the (001) interface with and with- 
out a Schottky barrier. The symmetry concepts still hold, 
since the least-decaying complex band of the tunneling 
barrier has the same symmetry properties as the conduc- 
tion band (both are of Ai character). We have seen that 
resonant interface states existing in the vicinity Ep for 
minority spin can be important for the current spin polar- 
ization. Even in the case of a Schottky barrier they can 
provide a resonant tunneling channel for the minority- 
spin electrons decreasing the spin injection effect. 

(iv) In the case of a higher bias, a wider energy range 
must be considered, so that states with k|| ^ as well 
as resonant states close enough to Ep become relevant. 
Then the spin polarization decreases, but spin injection 
is still achieved. 

It can be argued that a well ordered interface, neces- 
sary for our symmetry selection rule, is difficult to real- 
ize experimentally. However, recent successful attempts 
in increasing the spin injection efficiency have been re- 
ported, accompanied by an improvement of the quality 
of the interfaceii^ Therefore, we are optimistic that our 
work will motivate further research in this direction. 
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FIG. 16: 2DBZ for Fe(OOl) adapted to the Fe/SC(001) inter- 
face (thick line). For comparison the 2DBZ for a free Fe(OOl) 
surface (thin solid line) and the (001) cut through the bulk 
Brillouin zone (thin dashed line) are also shown. The vec- 
tor along the [001] direction varies between ±27r/cipc. 



APPENDIX A: GEOMETRY OF THE REDUCED 
2DBZ IN (001), (111), AND (110) 

First we investigate the backfolded bands for the 
Fe(OOl) half-space. In Fig. [H the 2DBZ and the cut 
perpendicular to the [001] direction of the bulk Brillouin 
zone are shown. At the F point two bands are available: 
the not backfolded band corresponding to the bulk band 
along F-H, i.e., the A direction. This direction has a 
C4t, symmetry. Also one band is backfolded by apply- 
ing once the reciprocal surface lattice vectors bi or b2, 
corresponding to the bulk band along N-P-N, i.e., the D 
direction. 

Next we discuss the (111) case. In Fig. [13 the (111) 
cut through the bulk Brillouin zone and the 2DBZ are 
shown. As for the (001) orientation additional bands are 
obtained by backfolding to the F point due to the large 
lattice constant in the Fe half-space. The not backfolded 
band at the F point corresponds to the bulk band along 
the F-P-H high symmetry line, because in the 2DBZ the 
kz vector extends from to \/3 in units of 2Tr/dpe. 

By applying bi or b2, together one additional band is 
backfolded to the F point. It corresponds to the bands of 
the same F-P-H high symmetry line as the not backfolded 
band. 

The 2DBZs of Fe and SC(llO) and the corresponding 
(110) cut through the bulk Brillouin zones are shown 
in Fig. ^1 In Fe(llO) there are four backfolded bands 
by applying the reciprocal surface lattice vectors bi and 
b2. The not backfolded band corresponds to the F-N 
bulk band (S direction). By applying the bi reciprocal 
lattice vector once, one ends up with a backfolded band 
that is not along a high symmetry line. More precisely, 
it is along the shortest line that connects two P points 
in the bulk Brillouin zone. Written in components of 
the reciprocal bulk lattice vectors, this backfolded band 



11 



ky[27C/C|re] 




corresponds to the bulk band along 





If both reciprocal lattice vectors, bi and b2, are applied, 
the backfolded band is parallel to the one, if only bi 
is applied, but shifted in the (110) direction. Written 
in components of the reciprocal bulk lattice vectors, it 
corresponds to 



FIG. 17: (111) cut through the bulk Brillouin zone (dashed 
line) and 2DBZ (sohd hnc) of Fc(lll). The kz vector varies 
in the 2DBZ between ±\/3 x 2n/dFe- 
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FIG. 18: 2DBZ (solid lines) and the (110) cut through the 

bulk Brillouin zone (dashed lines) for Fe(llO) (left panel) and 
SC(llO) (right panel). Values are given in 2n/d, with d being 
dpe and dsc for the left and right panel, respectively. 





The next backfolded band is obtained by applying b2 
once that corresponds to the bulk band between the high 
symmetry points N-H. In summary due to the large two- 
dimensional surface unit cell and small 2DBZ there are 
four bands available at the F point for the (110) orienta- 
tion in the Fe half-space including two bands correspond- 
ing not to a bulk band along a high-symmetry line. For 
these bands the only symmetry operation is the identity 
one, so that they are allowed to couple to any SC bands 
due to symmetry. 

In the SC (110) SBZ the not backfolded band corre- 
sponds to the bidk high symmetry line T-K (E-dircction). 
By applying bi one additional band is backfolded to the 
r point, exactly on the not backfolded band. So the 
maximum conductance can be 2 e'^/h per spin channel. 
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